On-site correlations in optical lattices: band mixing to coupled quantum Hall puddles 
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We extend the standard Bose-Hubbard model to capture arbitrarily strong on-site correlations. In 
addition to being important for quantitatively modeling experiments, for example, with Rubidium 
atoms, these correlations must be included to describe more exotic situations. Two such examples 
are when the interactions are made large via a Feshbach resonance, or when each site rotates rapidly, 
making a coupled array of quantum Hall puddles. Remarkably, even the mean field approximation 
to our model includes all on-site correlations. We describe how these on-site correlations mani- 
fest themselves in the system's global properties: modifying the phase diagram and depleting the 
condensate. 



Optical lattice systems, where a dilute atomic gas is 
trapped in a periodic potential formed by interfering laser 
beams, provide a close connection between solid state 
systems and atomic physics f\\. The models used to de- 
scribe these systems generally assume that each lattice 
site's wavefunction is easily built up from single parti- 
cle states [2]. Here we argue that this approximation 
is inappropriate to quantitatively model current experi- 
ments, and sometimes fails more drastically, e.g., for res- 
onant bosons. We show how to include arbitrary on-site 
correlations via a generalized Hubbard model, which can 
be approached by standard methods. By construction, 
the mean field approximation to our model captures all 
on-site correlations. 

Our method's key idea is to first consider deep lat- 
tices, where lattice sites are isolated and then solve the 
few-body problem on each site. Next, truncating to this 
few-body problem's low energy manifold, we calculate 
how tunneling couples the few-body states on neighbor- 
ing sites. The resulting theory resembles a Hubbard 
model, but with number-dependent hopping and inter- 
action parameters. We show that the corrections to the 
ordinary Bose-Hubbard model captured by this theory 
are crucial to quantitatively describe current Rubidium 
experiments. They become even more important when 
the 3D scattering length a becomes a significant fraction 
of the size of the Wannier states such as in recent exper- 
iments on Cesium atoms near a Feshbach resonance [3|. 
This approach is also essential to describe more exotic 
on-site correlations; as one example, one can rotate each 
lattice site, creating a lattice of coupled "quantum hall 
puddles" [1, m . Related ideas can be applied to double 
well lattices and coupled "plaquettes" of four sites 0] • 
We explore the impact of the on-site physics on the ex- 
tended system's phase diagram. 

Our approach is most simply illustrated by a single- 
component Bose gas in a cubic sinusoidal lattice potential 
Vp{x, y, z) = Vq ^r]=x,y,z sin^(7r7//d) with Hamiltonian 



where m is the particle mass, /i is the chemical potential, 
and ip and ip^ are bosonic annihilation and creation oper- 
ators. Adding an additional trapping potential presents 
no additional difficulties. 

Constructing the effective Hamiltonian. — For each iso- 
lated site, we proceed to build up the many-body states 
from the solution of the n-body problem at site j: 
(rir-- ,rn||n)j_= V^n(ri-Rj,--- , rn -Rj), which obeys 
— ^n|^)j where Hj is the same as Eq. ([T])'s i^/, 
except replacing the periodic potential Vp there with an 
on-site potential Vj. A convenient approximation is to 
take Vj{r) = {mu;^ /2){r - Rj)^ with u; = 2^/VoE^/h, 
the harmonic approximation to the site located at Rj. 
For each site filling n, we restrict our on-site basis to the 
lowest energy n-body state; however, including a finite 
number of excited states is straightforward. Note that 
even in the non-interacting case these states are not Wan- 
nier states. The principle difference is that states defined 
in this way are non-orthogonal. From these, however, 
one can construct a new set of orthogonal states |n)j, 
which hold similar physical meaning. In the noninteract- 
ing limit, the \n)j approximate the Wannier states. 

Because the single-site wavefunctions decay like Gaus- 
sians, it typically suffices to build up the effective 
Hamiltonian from neighboring sites. In particular, con- 
sider two sites L and i?, and the space spanned by 
Wl^tir) = \nL)L ® \nR)R, with overlaps = 
|m + l,n — 1). To lowest order in the over- 
we can define orthogonal It il^tir) by taki ng 
{1 /2) [^^^^'^^^ 
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Within this restricted basis, the effective 
Hamiltonian for these two sites is i^eff = 
T.n,m,n',m' W ,m') {u' , TTi' \H f\n, m) (n,m|. Evalua- 
tion to lowest order in S'^^^^ yields on-site energy terms 
5^n,m(^n Em)\n,m) {n,m\ and a "hopping" term 
-X^^^t(^^)|m + l,n- l)(m,n| +H.C. with 
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En = {n\Hf\n), 
t^^^) = -(m + l,n- l|i^/|m,n) 
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Additionally there is an interaction term U 
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to 0(5^), consistent with the rest of our calculations. In 
the remainder of this paper we will neglect the off-site 
interaction, Eq. (|3]), and the last term in Eq. ([2]). The 
former is rigorously justified as it falls of exponentially 
faster than the other interaction terms. Formally, the 
non-orthogonality contribution to Eq. ([2j) is suppressed 
only by a factor of (Vo/ERy^'^ with Er = h^iT'^/{2m(P), 
but as shown in Fig. [TJ it is typically small. 

The simplest many-site Hamiltonian which reduces to 
this one in the limit of two sites is 

H = - i™"^|m+l)Jn-l)^.H,H. 

;m,n 

+ ^i?„|n),(n|,, (4) 

where ^i^ij) indicates a sum over nearest neighbors i 
and j. At higher order, one generates more terms such 
as next nearest neighbor hoppings, pair hoppings, and 
longer range interactions. 

Calculating the Hamiltonian parameters. — Here we 
consider the cases of weak interactions, resonant inter- 
actions, and coupled quantum Hall puddles. 

In the limit of weak interactions, one can esti- 
mate the parameters in Eq. (|4j) by taking the on- 
site wavefunction to be ifju oc exp(— ^J^-,^ r|/2cr^), 
with variational width To leading order in 

a/d we find E^ = Er[(?>^V^/E^ - fi/ER)n + 

Unin-1) (l - ^{a/d){n - l)(Fo/^i?)'/^) ] with U = 

{a/d)V2^{Vo/ERf^^ and t^^^) = t^^^^- 

^%^(m + n - l){y^/ERfl^] with 



t = Fo(^V4 - 
l)g-(^V4)VVo/£;K^ Note that as expected, interaction 
spreads out the Wannier functions, increasing the t'^^^)'^ 
and decreasing the E^s. Fig. [T] shows several of the re- 
sulting t^^^^ as a function of Vq for parameters in typi- 
cal optical lattice experiments with ^^Rb. Also shown 
are t*^^^^ and the next-nearest-neighbor hopping, tnnn, 
calculated from the exact Wannier states. Our esti- 
mates are consistent with previous work regarding fsn- 
dependence [1, 0, [13] , validating our approach. As can be 
seen in Fig.fT^c), the relative size of the nearest neighbor 
hopping tnnn/t is 10% (1%) for Vo = 3Er {Vo = IOEr), 
justifying our approximation of including only nearest 
neighbor overlaps to describe the system near the Mott 
state. We also see that even for this weakly interacting 
case the number dependence oft is crucial for a quantita- 
tive description of the experiments. Similarly, the num- 
ber dependence of the on-site interaction is quantitatively 
significant. This latter deficiency of the standard Hub- 
bard model has been noted in the past, for example by 
the MIT experimental group [Tlj . 
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FIG. 1: (color online) (a) On-site energy with non- interacting 
energies subtracted off, In = En — {3h(jj/2 — /i)n, and energies 
scaled by Er — H^tt^ /{2md^), using typical ^^Rb parameters: 
lattice spacing d = 532nm, scattering length a = 5.32nm. 
Dashed: neglecting on-site correlations, solid: including cor- 
relations, (b) Representative hopping matrix elements with 
on-site correlations, relativ e to those neglecting on-site corre- 
lations, r^^^^ = t(™)/t^(m + l)n, as a function of lattice 
depth Vo on a scale. Bottom to top line: t^^^^ , t^^^^ , t^^^^ , t^^^^ . 
(c) For comparison, t'^^^^ /Er calculated from the exact Wan- 
nier states (upper curve) along with our Gaussian approxima- 
tion to it with and without non-orthogonality corrections (sec- 
ond and third highest, respectively); also shown is the next- 
nearest neighbor hopping matrix element (bottom curve). 
The effective Hamiltonian parameters are calculated pertur- 
batively 'm. a/d for a Gaussian ansatz. 



For more general experimental systems, one needs to 
include still more on-site correlations. As our first ex- 
ample, we consider lattice bosons near a Feshbach reso- 
I2I . 13], describing, for example, ongoing Cesium 



nance ^ 

atom experiments [3[. We restrict ourselves to site oc- 
cupations n = 0,1,2, for which we have exact analytic 
solutions to the on-site problem for arbitrary a in terms 
of confiuent hypergeometric functions [13] • Fig. [2] shows 
graphs of E^ and rescaled by huj as a function of a 

in the deep lattice limit. 

As Fig. [2](b) illustrates dramatically, the hopping from 
and to doubly-occupied sites is strongly suppressed near 
the Feshbach resonance when atoms occupy the lowest 
branch, and is enhanced for the next-lowest branch. The 
former has iniplications for studies of boson pairing on 
a lattice [12|, ll3[, showing that one must dramatically 
modify previous models near resonance, and, as will be 
discussed more below, the latter implies a substantial re- 
duction of the n = 2 Mott lobe's size for repulsive bosons. 

We give one further example, namely the case, sim- 
ilar to the one discussed in 0, [HI, where the indi- 
vidual sites of the optical lattice are elliptically de- 
formed and rotated about their center. This is ac- 
complished by rapidly modulating the phase of the 
optical lattice lasers to generate an appropriate time- 
averaged optical potential. At an appropriate rota- 
tion speed Vt the lowest energy n-particle state on 
each site is a = 2 Laughlin state ?/;n(ri, . . . , Tn) = 

[nr<,=iK -^J-)')] e-^.(l-^l'+^?)/(^^'\ where we 
define Wj = Xj + iyj, and ATn is a normalization fac- 
tor with phase chosen to gauge away phase factors ap- 
pearing in t^^^) , Truncating to this set of states for 
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FIG. 2: (color online) Left: On-site two-particle energy as 
a function of scattering length a rescaled by the on-site har- 
monic oscillator energy hoo = 2\^Vo Er, for the two lowest 
energy branches. The corresponding characteristic length is 
i — yjh/muj. Right: Log plot of rescaled hopping matrix 
elements t(™)/(^(m + l)nt). Solid and dashed curves are 
t(^^V(^^) and t(^^V(2t), respectively. We have chosen the 
lattice depth Vb = 15^i?; this affects only the horizontal scale. 
In the ordinary Bose-Hubbard model, t*^™^/( ./(m + l)nt) = 
1 for all m, n, as confirmed by this figure's a — 0~ lowest 
branch and a = 0^ second branch limits. The resonance at 
—d/a — separates the molecular side (left), from the atomic 
side (right). 



n = 0,1,2, we produce an effective Hubbard model 
of the same form as Eq. (|4]). The hopping parame- 
ters for asymptotically deep lattices Vq/Eji ^ 1 are 
^(01) _ ^(02) _ t{7ry32){Vo/ER)^/^, and t^^^) = 
t{7r^/1024){Vo/ Er) where Er = h^n^ / {2m(P) is the re- 
coil energy, and t is the same as in the weakly interact- 
ing case treated above; the interaction parameters are 
Era = ^^^^^m{m — 1) — jJ^m. One particularly interesting 
aspect of this model of coupled quantum Hall puddles is 
that when the system is superfluid, the order parameter 
is exactly the quantity defined by Girvin and MacDon- 
ald [il, llol to describe the nonlocal order of a fractional 
quantum Hall state. Thus when one probes the super- 
fluid phase stiffness, one directly couples to this quantity. 

Mean-field theory. — The true strength of our approach 
is that the resulting generalized Hubbard model is 
amenable to all of the analysis used to study the standard 
Bose-Hubbard model. In particular, we can gain insight 
from a Gutzwiller mean-field theory (GMFT) 0, [l^. 
This approximation to the ordinary Bose-Hubbard model 
gives moderate quantitative agreement with more so- 
phisticated methods: for example, the unity site filling 
MI/SF transition on a 3D cubic lattice occurs at {t/U)c = 
0.03408(2) while GMFT yields (t/U)c = 0.029 (l8|. 

In the ground state |^), we introduce mean fields 
^rn = (^|m + 1) (m|^). Neglecting terms which are 
quadratic in SL]^ = |i,m + l)(i,m| 
tonian is Hmf = Hmf^i with 



^rn, the Hamil- 



Hmf.i En, \n). {n\, 




Cm \m - l)i (m| ■ - Cm'fm-l + H.C. 



, (5) 



FIG. 3: (color online) Representative Gutzwiller mean- field 
theory phase diagrams, showing constant density (black, 
roughly horizontal) and constant ^ = C1+C2 + C3 (red, roughly 
vertical) contours: ^, similar to the condensate density, is a 
combination of the mean fields Cm, defined after Eq. (|5]). Den- 
sity contours are n = {0.01,0.2,0.5,0.8,0.99,1.01,1.2,...} 
and order parameter contours are ^ = {0.2,0.4,...}, ex- 
cept (d) where we take contours ^ = {0.02, 0.04, . . .}. The 
phase diagrams are functions of /ieff = /j/Erec and teff = 
exp where the lattice depth Vb is the natural 

experimental control parameter. We plot versus teff, instead 
of Vb, as this is closer to the Hamiltonian matrix elements 
and more analogous to traditional visualizations of the Bose- 
Hubbard phase diagram, (a) Ordinary Bose-Hubbard model 
for a = 0.01c/, (b) resonant lattice bosons restricted to fill- 
ings n = 0, 1, 2 with a = O.Old, on the next to lowest energy 
branch on the a > side of resonance, (c) resonant lattice 
boson model with a — d^ and (d) FQH puddle array model 
taking uo — = O.l^rec (see text for details). 



where H.c. denotes Hermitian conjugate, z is the lattice 
coordination number, and = Xlm • 

Truncating the number of atoms on a site to n < nmax, 
we self-consistently solve Eq. ([5]) by an iterative method. 
We start with trial mean- fields, calculate the lowest en- 
ergy eigenvector of the (n^iax + 1) x (nmax + 1) niean field 
Hamiltonian matrix, then update the mean-fields. We 
find that it typically suffices to take nmax roughly three 
times the mean occupation of the sites. Fig. [3] illustrates 
how the density dependence of the parameters introduced 
by the on-site correlations modify the GMFT phase di- 
agram — particularly the phase boundary's shape, and 
the density and order parameter in the superfluid phase. 

As one would expect, the topology of the MI/SF phase 
boundaries are similar to that of the standard Bose- 
Hubbard model, but the Mott lobes' shapes can be signif- 
icantly distorted. Within mean-field theory the bound- 
ary's shape can be determined analytically by taking 
1^) = |n - 1) + Vl - e'^ - e^/n-i \n) + e | n + 1), and 
expanding {^\Hmf\^) to quadratic order in e and e^ 
The Mott boundary corresponds to when the energy ex- 
pectation value's Hessian changes sign; this boundary oc- 
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FIG. 4: Single lobe of the Mott insulator/superfluid bound- 
ary. Complete characterization of the mean-field lobe shape 
in the zt-Jl plane for all possible t^'s (see Eq. (|7| for defini- 
tions), (a) Fix t~ = 0.5, vary from 1 (outer curve) to 21 
(inner curve) in steps of 2. (b) Fix = 1.5, vary t~ from 
(outer curve) to 1 (inner curve) in steps of 0.15. 

curs when 

The five scaled parameters 

En — En-1 En-\-l + En-1 — 2En 

-^{n,n) ^(n,n+l) ^(n— l,n) 

completely characterize the shape of a filling-n Mott 
phase boundary. Varying t and ft while fixing the other 
parameters then maps out a Mott-lobe like feature in the 
t and fl plane, is illustrated in Fig. [H 

Summary and discussion. — We have demonstrated a 
novel approach to strongly correlated lattice boson prob- 
lems. One constructs a model by truncating the on- 
site Hilbert space to a single state for each site filling n 
and includes only nearest-site, single particle hoppings. 
While this approximation captures many multi-band ef- 
fects, it is not a multi-band Hubbard model and in partic- 
ular retains the ordinary Bose-Hubbard model's simplic- 
ity. In this method, arbitrary on-site correlations may 
be treated, even in the mean field theory, and conse- 
quently it captures the condensate depletion, modified 
excitation spectra, altered condensate wavefunction, and 
altered equation of state characteristic of strongly inter- 
acting bosons. We have calculated the mean field Mott 
insulator/superfluid phase boundary analytically and ob- 
servables across the phase diagram numerically. 



Finally, although we have truncated to a single many- 
body state for each filling n, no difficulty arises from 
including on-site many-body excitations in the Hamilto- 
nian. These are especially important, for example, for 
double well lattices and spinor bosons. The ideas also 
extend straightforwardly to fermions; see Refs. [l^, 
[21I . [22! for related considerations. 
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